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Abstract 

This review is devoted to collecting some results on the high spin expansion of (minimal) anomalous 
dimension. Thanks to the recent rationale on integrability, planar N = 4 super Yang-Mills theory 
(or its AdSs x S 5 string counterpart) represents a very practicable field. Here the attention will be 
restricted to its sector of twist operators, although the analysis tools are quite general (in integrable 
theories). Some structures and ideas turn out to be general also for other sectors or gauge theories. 
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1 Framework and beyond 



We will move our investigation in within the maximally supersymmetric gauge theory in planar limit, 
i.e. for number of colours iV — > oo and coupling g YM — > 0, so that the 't Hooft coupling 



may stay finite. Among the different sectors (perturbatively closed under renormalisation), we also pick 
up the twist sl(2) sector, spanned by local composite operators of trace form 



where T> is the (light-cone) covariant derivative acting in all the possible ways on the L complex bosonic 
fields Z. Trace ensures, of course, gauge invariance. The Lorentz spin of these operators is s and L is the 
.R-charge which also coincides with the twist (classical dimension minus the spin). Besides, this sector 
may be described - thanks to the AdS/CFT correspondence pQ - by spinning folded closed strings on 
AdS 5 x S 5 spacetime with AdS 5 and S 5 angular momenta s and L, respectively [21 E]- 

As being in a conformal model, suitable superpositions of operators form dilatation operator eigen- 
vectors with definite dimensions (eigenvalues), which are made up of a classical part plus an anomalous 
one. For instance, in the sector ( 11. 21) this spectral problem shows up dimensions 



where , y(g,s,L) is the anomalous part. According to the AdS/CFT strong/weak coupling duality, the 
set of anomalous dimensions of composite operators in M = 4 SYM coincides with the energy spectrum 
of the AdS 5 x S 5 string theory ([H [21 E| and references therein), although the perturbative regimes 
are interchanged. The highly nontrivial problem of evaluating the anomalous part in M = 4 SYM 
was greatly simplified by the discovery of integrability in the purely bosonic so(6) sector at one loop 
[I]. Later on, this fact has been extended to all the gauge theory sectors and at all loops in a way 
which shows up integrability in a weaker sense, but still furnishes the investigators many powerful tools 
[5]. More in detail, any operator (e.g. of the form (11.21) ) has been thought of as a state of a 'spin 
chain', whose hamiltonian is, of course, the dilatation operator itself, although the latter does not have 
an explicit expression of the spin chain form, but for the first few loops. Nevertheless, the large size 
(asymptotic) spectrum has turned out to be exactly described by certain Asymptotic Bethe Ansatz-like 
equations (the so-called Beisert-Staudacher equations, cf. [5JE] and references therein). In other words, 
the anomalous dimensions coincide with the energies given by the Bethe Ansatz solutions (or roots): 
this is, of course, a great simplification of the initial spectral problem. 

Unfortunately, this works only for infinitely long operators: anomalous dimensions of operators with 
finite length depend not only on Asymptotic Bethe Ansatz (ABA) data but also on finite size 'wrapping' 
corrections: in the perturbative expansion wrapping effects are observed [7] starting from the order g 2L . 
Recent progress has shown that a set of Thermodynamic Bethe Ansatz (TBA) equations [8] or, together 
with certain additional information, an equivalent F-system of functional equations [9] provides a basis 
for exact (any length at any coupling) predictions for anomalous dimensions of planar M = 4 SYM. 



A = g 2 YM N = SvrV 






A (g, s,L) = L + s + j(g, s, L) , 
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However, in the sl(2) sector of = 4 SYM relevance of wrapping effects seems to be reduced, even 
for short operators, if one goes to the high spin limit. For instance, findings of [10] have showed that, at 
least at twist two and up to five loops, wrapping corrections start contributing at order O ((In s) 2 /s 2 ). 
This fact has pushed the idea of applying ABA techniques to the study of the high spin limit of twist 
operators, which - on the other hand - had already received much attention in the past literature. 
Indeed, the high spin behaviour of the anomalous dimension 

l(g,s,L) = A(g,s,L)-L-s, (1.4) 

shows a Sudakov behaviour 

j(g,s,L) = f(g)\ns + ... , (1.5) 

determined by the so-called universal (since it does not depend on L or the flavour) scaling function, 
f(g) [Til [12| [T3| [61 [H] . Actually, this behaviour is more general than in planar M = 4 SYM and the 
adjective scaling is due to the linear value of the cusp anomalous dimension with coefficient f(g)/2 while 
the cusp angle tends to infinity [15] 0. This large angle behaviour is due to the dominance by the lowest 
twist (= 2) in the renormalisation of the vacuum expectation value of a cusped Wilson loop with a very 
large angle. In the end, for an infinite angle cusp (i.e. with one light-cone segment) f(g) equals twice 
the cusp anomalous dimension of a light-cone Wilson loop [17]. Additionally, the one- loop problem 
(and thus /(<?)) stays exactly the same for twist operators in QCD as long as the partonic helicities are 
aligned [HI H5] ; this fact also has justified partially the great interest on the twist operators (11.21) . 

In general, the high spin limit of anomalous dimensions of twist L operators goes on as a series of 
logarithmic (inverse) powers^], 

oo 

j(g, s, L) = f(g) In a + f si (g, L) + £ 7 (n) (<?, L) (In s)~ n + O ((In s)" 00 ) , (1.6) 

n=l 

i.e. looks, at this order, like an expansion in the large 'size' Ins. Recently [6], the leading term (i.e. 
f(g)) was obtained - in the hypothesis of being wrapping free - from the solution of a linear integral 
equation directly derived from the ABA via the root density approach. Moreover, f(g) was carefully 
studied and tested both in the weak [121 E] and strong coupling limit [201 12H 1221 [23] . 

The sub-leading (constant) contribution f s [(g,L) received also much attention. In the ABA frame- 
work, it was shown [21] to come from the solution of a non-linear integral equation (NLIE). Then, it 
was obtained [25] starting from a linear integral equation (LIE). Explicit weak and strong coupling 
expansions performed using the LIE of [25] are present in [261 12Z] and agree with some string theory 
computations |28^l . Importantly, after results of [10] (cf. above) it may be inferred that both f(g) and 
f s i(g, L) are exactly given by this approach based on the ABA (without wrapping corrections). Besides, 
both terms fix the 1/s coefficients via the reciprocity relation which may still be, consistently with [TU] , 
wrapping free. 

1 Polyakov noticed as first that for cusped Wilson loop vacuum expectation value the charge renormalisation is not 
enough as in the non-cusped case, because of an extra logarithmic divergence due to the high bremsstrahlung at the cusp. 
He was led to consider cusps because of their importance in the loop dynamics (in euclidian space-time) [16] . 

2 With O ((In s)~°°) we indicate terms going to zero faster than any inverse powers of In s. 

3 Agreement between ABA [29] and string theory [28] results is now established up to two loops. 
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The latter reasoning would support that even sub- logarithmic terms from ABA be exact, but it is 
even more unclear if and at what extent this might be trusted. In [30] we studied and computed in a 
systematic way the ABA contribution to them, by using a set of integral equations. Possible wrapping 
corrections to our results are still to be determined. In this respect, in addition to TBA findings, also 
results on the string side of the correspondence in the spirit of [28] could be of fundamental importance. 
In fact, for instance, the one-loop contribution to the 1/lns term in the second of [28J comes partially 
from the AdS§ modes (wrapping free) and partially from the S 5 ones (wrappingjj. Nevertheless, a 
better understanding might come from exact (maybe TBA) calculations. 

Wrapping effects should be negligible^ if one takes the infinite twist limit 

L — 2 

s — > oo , L — »• oo , j = — fixed , (1.7) 

In s 

and restricts to the study of the scaling functions f n (g), fn\g) appearing in the expansion 

oo oo oo 

-y{g, s, L) = In 8 fn(9)f + E( ln s )~ r E fn\g)f + O ((In *)-") . (1.8) 

n=0 r=0 n=0 

For this reason, a reasonable amount of activity was also devoted to the study of limit (jl.7p . using ABA 
techniques. Results concerning f n (g) at weak coupling are present in [24]. Strong coupling behaviour of 
fn(g) was studied in [2H E21 E31 EH] by relying on linear integral equations. A detailed study of fn\g), 
r > 0, can be found in [35| 130] . 

In this review we want to give a summary of our activity in the framework of high spin limit of twist 
operators. We will first give a general description of the method we have used and whose main tools are 
integral equations derived from the ABA equations. Then, we will briefly report the most important 
results we obtained. 

This review is organised as follows. In Section 2 we review the ABA equations for twist operators. 
We illustrate the properties of the minimal anomalous dimension operators and explain a technique 
which allows exact computations of ABA contributions to the anomalous dimension by using the so- 
called Non-Linear Integral Equation (NLIE). In Section 3 we specialise to the high spin limit and show 
that, if one neglects O ((Ins) -00 ) terms, asymptotic anomalous dimension can be computed by relying 
on integral equations. In Section 4 we study the high spin limit at fixed twist. In Section 5 results in 
the scaling limit (jl.7j) are discussed. 



2 All-loop ABA and the (N)LIE 



Let us recall the ABA equations JSJ E] for the si (2) sector, 

^ ± i) L ( 1 + ^) = f[ e»*»*», (2.11 



4 We wish to thank N. Gromov who pointed out this double origin as he sees it clearly from the algebraic curve approach 
which reproduces one- loop ABA plus some 'vacuum oscillations' (wrapping). 

5 This should be true at least for small values of g. But, again, this reasoning is not rigorous. 
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where 



± 



u 



%i = x (u k ) = x(u k ± i/2) , x(u) = - 



u- 



X 



3 2 2 



A being the 't Hooft coupling. The so-called dressing factor [SHI EH IE] #(w, v) is given by 

oo oo 

r+l+2 i y(5 , )[gr(«)gr+l+2 i /(f) - (u)<?r+l+2i/ («)] , 



(2.2) 



(2.3) 



r=2 i/=0 



the functions /3 r)r+ i +2 i,(<7) = g 2r+2v 2 2 1 r u c rir+1+2 v(g) being 

°° 2r+2u+2fi I 



l){r + 2u) 



fl=U 



2/i 



2^+1 
— r — i/ + 1 

and q r (u) being the density of the r-th charge 

r-1 

s i I i \ 

q r (u 



2/i+ 1 
\i — v 



C(2/i + i; 



x (n) 



r-1 



(2.4) 



(2.5) 



It is now clear that configurations of Bethe roots, i.e. solutions of ( 12.11) . and the corresponding eigen- 
values of the energy are related respectively to composite operators and their anomalous dimensions in 
the sZ(2) sector of M = 4 SYM. 

In the si (2) sector states of twist L are described by an even number s of real Bethe roots which 
satisfy (12. ip . Bethe roots localize in an interval [—b, b] of the real line. In addition to Bethe roots, also 
L real 'holes' [TTj I3"5| El I2H 125] are present (a better explanation of the nature of holes will be given in 
the following). For any state, two holes reside outside the interval [—6, b] and the remaining L — 2 holes 
lie inside this interval. We will indicate with u£ these 'internal' holes. 

In this paper we will focus on the minimal anomalous dimension state. For such a state the positions 
of both roots and holes are symmetric with respect to the origin. For what concerns the internal holes, 
they all concentrate near the origin, with no roots lying in between. 

In general, an efficient way to treat states described by solutions of a non-linear set of Bethe Ansatz 
equations consists in writing a non-linear integral equation, which is completely equivalent to them. 
The non-linear integral equation is satisfied by the counting function Z(u), which in the case (12 .ip reads 
as 



Z(u) = $(u) 



s 

fc=i 



where 



$(u) 



[u,v] 



-2L arctan 2u — iL In 



[u,u k 



2x-(uf 



1 + 



2x+{uY 



2 arctanfn — v) 



2i 



In 



2x+ [u)x~ (v) 



1 - 



i9(u, v) 



2x~ (u)x+(v) 



(2.6) 

(2.7) 
(2.8) 
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It follows from its definition that the counting function Z( function of the real variable u, is a 

monotonously decreasing function and that 



lim Z{u) = tk(L + s) . (2.9) 

u— >±oo 

Therefore, there are L + s real points Vk such that e* 2 ^ = (— 1) L+1 . It is a simple consequence of 
the definition of Z(u) that s of them coincide with the Bethe roots. The remaining L points are called 
'holes' and their role will be of fundamental importance in what follows. As we anticipated before, 
for the minimal anomalous dimension state the internal holes concentrate near the origin, i.e. their 
positions u£ are determined by the relations 

Z(uf) = n(2h+ 1 -L), h = l,...,L-2. (2.10) 

Excited states are obtained by making different choices for the 'quantum numbers' h. It follows from 
the structure of (I2.10p that, for any state, depend in a non-linear way on the counting function 
Z{u). 

The non-linear integral equation for Z(u) is written by using a modification of the standard ideas 
underlying the procedure concerning the excited state NLIE [39], this modification being dictated by 
the physical situation with two 'important' external holes. Suppose that in the interval [—b, b] of the 
real line s Bethe roots and L — 2 holes are present. Then, by using Cauchy theorem we can express a 
sum over the Bethe roots of an observable 0(u) as 

yO(u k ) = - I" ^0(v)Z'(v) + lm I" —0(v-te)^-\n[l + (-l) L e iZ ^] + 
^ J-b 271 J-b 71 dv 

+ Im / " ^-0(-b + ty)^-\n[l + (-l) L e iZ ^ b+iy) }+ (2.11) 

+ Im f ^0{b + iy)^-\n[l + {-l) L e^ b+ ^]-Y J 0{uf). 
J-e 71 dy h=i 

The right hand side of (12.1 ip does not depend on e > as far as no poles of the integrands 0(w)-^- ln[l+ 
{ — \) L e tZ ^] lie in the region |Imit?| < e, |Rew;| < b. In addition, e must be kept sufficiently small, in 
such a way that 

|e iZ(z) |<l, (2.12) 

where z belongs to the integration contour of the last three integral terms in ( 12. lip . This condition is 
assured by the monotonicity of the counting function, i.e. in our case Z'(v) < 0, v 6 [—6, b]. We apply 
( 12. lip to the sum over the Bethe roots contained in (12.61) : 



Z{u) = —<f)(u,v)Z\v) + Y<f)(u,u < i ) )-lm[ —<f){u,v-ie)^-ln[l + (-l) L e iZ{v - le) }- 

J-b 2vr ^ J^ b 7T dv 

(2.13) 

l m f E ^{u, -b + iy)4- ln[l + (~1) L e iZ ^ b+iy) ] -Im [ ^-<j>{u, b + iy)^- ln[l + (-1) L e lZ{b+iy) ] . 
Jo 7T dy J_ e it dy 
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What we have obtained is a non-linear integral equation - for the counting function Z(u) - which 
describes - in a way which is alternative to the Bethe Ansatz equations - the minimal anomalous 
dimension state. Since Bethe roots are localised in an interval of the real axis (i.e. b < +00), NLIE 
(I2.13p is different from the equation introduced in [39], which contains integrations over the entire real 
axis. This difference reveals crucial in the specific problem of twist operators in the si (2) sector. Indeed, 
important simplifications to the structure of (I2.13P arise in the high spin limit: if we decide to neglect 
terms going to zero faster than any inverse power of In s, then it is possible to get rid of the last three 
non-linear terms in ( I2.13p . We are now going to show this important property. 

Coming back to (12. lip , we first use (12. 12ft in order to replace all the ln[l + (— l) L e iZ ( z >] with 



E(-i 

n=l 



n+l(_l) 

s 



. We obtain: 



dv 



k=l 



dv 



-OMZ-W + Im/ („ -«) £(-! 



n+l. 



-1) 



nL ( 



inZ(v—ie) 



+ Im/ ^Qi-b + iy)^- 1 



d 



n+l. 



dy 



■1) 



n=l 
D inZ(—b-\-iy) 



n 



+ 



+ 



n=l 



n 

fod+^B-n-i)" 15 - — E°(«i: 

e ^ " y n=l fc=l 



(2.14) 



In order to evaluate the non-linear terms in this expression, 

»6 



Im 



+ Im 



6 7T dV Z — ' 



n+l/ 



+ 



n=l 



n 

ginZ(— 6+ij/) 



7T 



+ 



71=1 



+ Im [°^-0(b + iy)^-f2(- r 

J-e 7T aj/ x 



71+1/ 



\7lL" 



n 

-,inZ (b+iy) 

n 



we first assume that it is possible to exchange the series with the integrals. Then, we use the following 
formula? 



dvOiv — ie 



d_ 

dv 



^inZ(v—ie) 



' i \ t d" 



E 

fc=0 
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" dyO(±b + iy)^e mZ{±b+iy) 
dy 



E 

k=0 



dx k 

k d k 



0[Z-\x) 



dx k 



x=Z(v—ie) 

0[Z-\x)] 



(2.15) 
(2.16) 



x=Z(±b+iy) 



We need to remark that results (12.151 12.16j) are correct if the above infinite sums make sense either as 
convergent or asymptotic series. In our case, we can use (12.151 12. 16j) . since the series we will get are 
asymptotic. 
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Using ( 12.151 12.16P in the evaluation of the above non-linear terms, the dependence on e cancels out, 
as it should be, and we are left with 



NL = J2 



<-i)» +i (-i) 



nL r 00 -2k 



n=l 



irn 



x=Z(b) 



d 2k i 2k d 2k+1 

fc=0 k=0 x=Z(-b) 

(2.17) 

Now, we are allowed to choose b in such a way that e lZ( - ±b ^ = (— 1) L : therefore, since 0(f) is bounded, 
all the terms in (12. 17ft proportional to the sine function are zero. Thus, we are left only with terms 
containing the cosine function, i.e., after summing over n, 



NL 



2 , m,. , 0M-D2/C+2 



fc=0 



(2fc + 2)!' 



2/ Ldx 2fc 



-0(^ x (x)) 



x=Z(6) 



x=Z(-6) 



(2.18) 



where B^{x) is the Bernoulli polynomial. Relation (I2.18P allows to write the final formula: 



fc=i 



L-2 



,2fc+l 



2vr 



h=l 



k=0 



(2k + 2)\ 



2/0+2 



5 



Q x 2k+1 



0(Z-\x)) 



x=Z(b) 
x=Z(-b) 

(2.19) 



The following expressions 



d 



x=Z(b) 

x=Z{~b) 
x=Z(b) 

x=Z(-b) 
x=Z(b) 

x=Z(-b) 



P'(b) 

[Z>{b)f 

0"'(b)-30"(b)^§-0'(b)^ + 30'(b)m 



{z'{b)y 



[Z'{b)\* 



give an idea of the form of the first nonlinear terms appearing in (12.191) . Next step is to apply (I2.18P 
to the non-linear integral terms (the last three ones) contained in (I2.13p . We have to replace 0(v) with 
—<f)(u,v). From the form of (I2.18P we realise that such non-linear terms are proportional to derivatives 



d n 
dv r 



n > 1 



(2.20) 



±i» 



When s — >oo,6=|(l + OQ)) and for such derivatives 



d n 
dv n ' 



lu,v] 



o 



±b 



b n+l 



(2.21) 



On the other hand, for the derivatives of the counting function we can give the estimates 

d n 



dv n 



Zv, 



O 



±1, 



(2.22) 



Putting all together, we conclude that when the spin s — > +00 the non-linear integral terms contained 
in (I2.13P are 0(l/b) = O ((lns)~°°). Therefore, if we decide to neglect terms O ((his) -00 ) and to focus 
only on corrections 0((ln s)~ n ), we are entitled not to consider all the non- linear integral terms in (I2.13p . 
In this case we are left with 

Z{u) = *(«) + f ^cj>(u, v)Z\v) + J2 u h) + O ((In s)—) , (2.23) 
•'~ b 71 h=i 

and this equation has to be solved together with condition (12.101) . which fixes the holes positions 
in terms of the counting function Ziu). It is important to remark that within our approximations 
non-linearity with respect to Z[u) enters only through (12.101) . 

Equations (12.231 12.101) are our starting point for the study of the high spin limit. They will be 
worked out in next section. 



3 Integral equations for the logarithmic terms 

Let us consider integral equation (12.231) satisfied by the counting function in the high spin limit. Passing 
to derivatives, we define 

a(u) = Z'(u) . (3.1) 



We have 



a(u) = &(u) 



L-2 



b dv d yrA d 

(p(u,v)(x(v) + 



2tt du 



h=l 



du 



(n,4 4) ) + 0((ln S )"°°) 



(3.2) 



We now decompose such equation in its one loop (with index 0) and higher than one loop (with index 
H) contributions. We set $(u) = $o( u ) + ®h(u), (fi(u,v) = <p (u,v) + <Ph(u,v), cr(u) = cr (u) + cr H (u), 
where 



$o( u ) — — 2Larctan2n, $#(u) = — zLln 



(f)o(u,v) = 2 arctan(w — v ) , <Ph(u, v) = — 2i 



1 + 



2x-(uy 



1 + 



2x+(uY 



hi 



2x + (u)x (v) 



+ i9(u, v) 



2x~(u)x+(v) 

and where, after neglecting quantities which are O ((lns) _oc ), <Jq{u) and (Th{u) satisfy 



a Q [u 



dv d 



d 



2tt du 



h=l 



/ \ , , , , f dv d , , . . . f dv d . , . 
a H (u) = ® H {u) + ——<Ph{u,v)(j h {v)+ I ——(p H (u,v)ao 
J_ b 2ir du J_ b 2ir du 



(«) + 



1 b ^ Mu > vMv) + ^du] [ 
~ b h=l 



M ,^ ) )-0oK4 ij )]+O((ln S )-° o ) . (3.6) 



(3.3) 
(3.4) 



(3.5) 



In (I3.5[ 13 .6 j) the notation stands for the one loop component of the position of the h-th internal 
hole. In order to solve the one loop equation ( 13. 5p . we consider (see 3.52 of [22]) a function <Jq\u), 
whose Fourier transform^ reads 

Function ( 13. 8 p satisfies the following important property [25] . 

f&0 P+OO 

duf{u)a {u) = duf{u)a [ s) {u) + ((Ins)" 00 ) , (3.9) 

bo J — oo 

which allows to extend to the whole real axis the integrations involving <Jq{u): a look at many loops 
equation (13. 6p shows that this is just the property we need in order to try to solve it. 

It follows from results in [12] and from numerical checks that we can extend to the entire real axis 
the integrations involving cjh(u): what we are missing are O ((Ins) -00 ) terms. Putting all these pieces 
together, the equation satisfied by <Jh(u) is 

»+oo 

a H (u) = & H (u 



- + 0O 



/-\-C>0 1 J ^ 2 J 

^M^v)a H (v) + J2^[Hu,4 ) )-M^^)] + 0((lnsr-) . (3.10) 
"°° h=i 

It is now convenient to pass to Fourier transforms. We have 

*H(k) = ^e-^[l-J (V2gk)}, (3.11) 

Mk,t) = -Szvr 2 ^^ [J2 r(-lY +1 J r (V2gk)J r (V2gt) Sg * [kt) + 

i i r= i 

OO OO 

+ sgn(t)^^c r , r+1+2 ,(( ? )(-l) r+ ^J r _ 1 (v / 2^A;)J r+2v (v / 2^) - 

- J r ^{V2gt)J r+2u (V2gk)) \ . (3.12) 



3 We adopt the following definition for the Fourier transform f(k) of a function f(u): 

-iku 



f(k) = / due-^f{u) . (3.7) 
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Therefore, the Fourier transform of (13.101) reads: 



sinh ^ sinhyj-oo H l r=1 



1 — sgn(fct) 
— e 2 



oo oo 



+ sgn(t)^^c rir+1+2! ,(^)(-l) r+ ^^(j r _ 1 ( v / 2^)J r+2 ^(v / 2^) - J r -i(V2gt)J r+2 u(V2gk) 



r=2 u=0 



L-2 



a H (t) + a { s \t) + 27r^e itu ^ 



i,=i 



-2tt 



Inserting in (13.131) the expression (13.81) for (t), introducing the 'magic' kernel K(t,t'), defined in [6] 

as 

2 



_W L-2 

e 2 



2 sinh ^ 



[e" ifeu " 



(0 



O ((lns)-°°) . (3.13) 



^(t,0 



oo oo 



.71=1 fc = l 1=0 

and restricting to > 0, we finally get the integral equation, 



(3.14) 



a H (k) 



[1- J (V2gk)]-g 2 



sinh | 



sinh | j 



dte-?K(V2gk,V2gt) 



L-2 



_ * L-2 

e 2 



e 2 cos(ts/A/2) x r/ x e 2 ^-^ m ^-^ 

(47rln2)^(t) + 2?r ^ - ; m ^ cosfiZ£ J + 2?r costw^ } + 



2 sinh 



If I 



9 sinh 



ft=l 



+ 7T \ cos ku^ — cos ku 

sinh | f-f L h 

^ /i=i 



(3.15) 



Since we are interested in the computation of the anomalous dimension, we go back to the general 
formula (I2.19P for the evaluation of observables, 



k=l 



Uk 



dv 
2^ 



L-2 



0(v)Z'(v)-J20(4 ] )-Yl 



(2tt 



,2/c+l 



h=l 



k=0 



(2k + 2)\ 



Boi. 



2k+2 







dx 2k ' 



z O{Z~\x)) 



If we specialise to the energy (anomalous dimension), we have 



0(v) = e(v) 



x + (v) x (v) 



Again, in the large spin limit, one has (with n > 1): 

d" 



dv n 



e(v) 



O 



b n+l ' 



x=Z(b) 
x=Z{-b) 

(3.16) 



(3.17) 



(3.18) 



together with the estimate ( I2.22p for the counting function. It follows that the non-linear terms con- 
tained in the expression (13. 16[) for the energy are 0(1/6) = 0(1/ s). Therefore, if we neglect O ((Ins) -00 ) 
terms, we are allowed to work with only the linear expression 

L-2 



l(9,s,L) = ^e(u k ) 

k=l 



— e(v)Z (^ 
b 27r 



^e( M i l) ) + 0((ln S )-) 



(3.19) 



h=l 
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which we find convenient to write in terms of the one loop and higher than one loop densities: 

rba . L-2 



l(g,s,L) = ^e(w fc ) 

k=l 



±L e (v)* H {v)- f° ^e{v)a Q {v)-Y,<uf) + 0{{\ns)-) . (3.20) 
^~ &0 h=i 



Extending the domains of integration to the real axis does not give problems, since we are neglecting 
O ((Ins) -00 ) terms: we get 



l{g,s,L) = ^2e(u k ] 

k=l 



+00 



dv 
2tt~ 



e(v)a H (v) 



+00 



dv 
2^ 



L-2 



-5>(4V0((M-°°) • (3-21) 



h=l 



Passing to Fourier transforms we have 



e 2 



e(k) = 2gV2ir——J 1 (V2gk) 



(3.22) 



and, consequently, 



l(9, s,L) 



dk 



_ JA1 
e 2 



L-2 



, 4ir 2 a " k 
Comparing ( 13. 23ft with (13.130 . we see that 



2gV2ir——J 1 (V2gk) a H (k)+a ( s) (k)+2n cos kuf +0((lns)"°°) . (3.23) 



h=i 



y(g, s, L) = - Urn a H (k) + O ((In s)" 00 ) 

7T fc->0 



(3.24) 



which extends the Kotikov-Lipatov relation [40] to all the sublogarithmic 0((lns)~ n ), n > 1, con- 
tributions and allows to compute the high spin anomalous dimension from the higher than one loop 
density. 

For computational reasons, it is more convenient to use the function 



S(k) 



sinh 



\k\ 



TT\k\ 



a H (k) - tv 



_!M L-2 

e 2 



sinh M h 



cos ku^ — cos fcw^ 



7(3, s,L) = 21imS'(fc) .(3.25) 

A;->0 



The function (13.251) satisfies the integral equation (for k > 0) 

S(k) = -[1 - J (V2^)] - g 2 / -e-li^(v / 2^, v^) 
fc Jo 71 



1Z_ 
— P 2 

i ——rS{t) - 4tt In 2 6{t) - tt( L - 2) — 

I sinh 7T si 



sinh I 



2tt 



1 -i is 

1 — e 2 cos ^= 



sinh I 



+ 



(3.26) 



1 L-2 

e2 ^-^ r (-^ 

+ 7T > COStMt — 1 



^nh I 



l 7i 



Ag 2 lnsK{V2gk,0) + 4g / 



dt 



e* - 1 



K*(V2gk,V2gt) + 



+ -[1- J (V2gk)]+Ag 2 lE K(V2gk,0) + g 2 (L-2) / dte~^K(V2gk, V2gt) — 



1—62 



sinh I 



+00 



L-2 



- g 2 dtk(V2gk,y/2gt) 



h=l 



cos tu 



in 



sinh I 



dte~^k(V2gk, V2gt)^—rS(t) , 

sinh 2 



12 



where K*(t, t') = K(t, if) — K(t, 0). We remark that such equation depends on the position of the holes 
(but not on both , uff) and that the quantity is also an unknown of the problem and has to 
be determined by solving equation (12.101) . 

Z[uf) = 7r(2h + 1 - L) , /i = l,...,L-2, 

which introduces non-linear effects in the equation for the density S(k). Equation f!3.26|) is exact if we 
neglect - in the high spin limit - terms which are O ((Ins) -00 ) and is the main integral equation of this 
papei0. It will be our starting point in order to investigate the high spin limit of (minimal) anomalous 
dimension of twist operators in the s/(2) sector. 

We will study two cases. First, we will consider the limit: 

s — > oo , L fixed. (3.27) 

Then, we will focus on the scaling limit (11. 7\i [TT] : 

L-2 n , 

s — )• oo , L — > oo , j = — fixed . 

In s 

The case ( 13. 27ft will be reported in Section 4, the case (11. 7p in Section 5. 



4 High spin at fixed twist 

In this section we report our results on the fixed twist case (13.271) . We naturally make a move from 
equation ( I3.26p . in which the various contributions to the known ('forcing') term are separated according 
to their power of In s. As a consequence of the structure of the forcing term, the high spin expansion 
of the function S(k) goes on as a series of logarithmic (inverse) powers ( 11.61) : 

oo 

S(k) = S(n) ( k ) (M" n + O ((Ins) -00 ) . (4.1) 

n=-l 

Consequently, the anomalous dimension at high spin follows the same fate: 

oo 

j(g, s, L) = f(g) In s + f sl (g, L) + ^ 7 W G?, L) (In s)~ n + O ((In s)" 00 ) . 

n=l 

Coming back to S(k), the term proportional to Ins satisfies the BES linear integral equation, 

/+oo j. 
dte-^K(V2gk, ^9 t )—^T S[ ~ 1) ( t ) > ( 42 ) 

7 Actually, for twist two there is a simple way to get also the O (^p) and O Q) terms from the first of (|3.26|1 : it is 
sufficient to replace in the argument of the cos function the quantity , which represent the leading contribution to the 
position of the external hole at large s, with the more accurate estimate [JT] (l + -\- 0(1/ s 2 )) . We then get the 
same result as (29) of [2l)] . 
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whose solution determines the universal scaling function through f(g) = 2S'( _1 )(0). The BES equation 
was introduced in [121 E] and thoroughly studied in the weak [121 Ej and the strong coupling [201 
limit. The interested reader can refer to these papers for a detailed study of equation (14. 2p 



The four subsequent terms - independent of s - appear in the linear integral equation for the density 
which determines the virtual scaling function f s i(g,L). This equation was written in [25] (equation 
4.11), then it was re-obtained in [26] and used there and in our contemporaneous paper [27] (where it 
appears as equation 3.3). In notations used in this paper it reads (k > 0) 

S {0) (k) = V / -P-K*(y/2gk, yftgt) + 
Jo e — i 

+ j[l-J (V2gk)}+4g 2 lE k(V2gk 1 0)+g 2 (L-2)£ dte^ K(V2gk, ^gt)^^- - (4.3) 

- g 2 j^ dte^k(V2gk,V2gt)— ^°\t) , 

the virtual scaling function f s i(g,L) being 

f al (g 1 L)=2S®(0). (4.4) 

As in the case of f(g), an explicit expression for f s i(g,L), interpolating from weak to strong coupling 
has not been found yet. What we did is expanding equation ( 14. 3 p in a systematic way for small g, thus 
getting the weak coupling convergent series for f si (g,L) (formula (4.1) of [27]): 



f sl (g, L) = ( lE - (L - 2) In 2)f(g) + 8(2L - 7)C(3) (jj=\ + 
-| (A(3)(L - 4) + 3(211 - 62)C(5)) (j=^ + 

+ ^ (vr 4 C(3)(3L - 13) + 75(46L - 127)C(7) + 5(1LL - 32)vr 2 C(5)) f + 
1 28 

— vr 6 C(3)(llL - 49) + 8 (2695C(9)L + 16C(3) 3 L - 7156C(9) + 56C(3) 3 ) + 

+ f(25L - 64)^(7) + A(i03L - 310) A(5) J " + 

+ ^tt 4 C(7)(295L - 772) + ^tt 2 (1519C(9)L + 24C(3) 3 L - 3628C(9) - 88C(3) 3 ) + 
+8 (33285C(11)L + 536C(3) 2 C(5)L - 86082C(11) - 1728C(3) 2 C(5)) + 

^(2023L - 6266) A(5) + ^6L - 13231)A(3)\ / 9 \" + _ _ _ _ { ^ 



945 v ' ^ v ' 14175 J \^2J 

In addition, in [27] we performed the strong coupling analysis, by means of analytical and numerical 
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computations. For the first leading terms (formula (4.12) of [27]) we provided the following outcome: 



f al (g,L) = 2V2i 



, 2^2 3 In 2 , 2^2 c + (2-L)7T 
In ci 7= — m 



9 



2^2' 



9 



2V2ng 



K , 2^2 
m 



k-i 



8ir 2 g 2 



9 



2y/2i 



where C\ — 1, Co = 61n2 — ir, K = 0(2) is the Catalan's constant and (see also [26] 

4K-9(ln2) 2 



k-i 



AV2tv 2 



-0.0118253... . 



9 3 ' 
(4.6) 



(4.7) 



Importantly, in this asymptotic expansion we could show that the only trace of the twist comes up in 
the piece cq + (2 — L)n and thus cancels out completely, at order 0(s°), in the asymptotic (large g) 
expansion of A — s = 7 + L. It follows that the constant term (i.e. 0(g°)) in A — s = 7 + L at order 
0(s°) is 6ln2+n for any twist: this allowed successful comparisons with string theory results [28], which 
does not distinguish between null and small values of L. 

An alternative analysis of the strong coupling limit can be found in [26]. This computation was 
done adapting the method introduced in [22] for f(g). 

For what concerns j^ n '(g, L), i.e. the 0((lns) _n ), n > 1, contributions to the anomalous dimensions, 
they are 'driven' by the holes depending parts of (13. 26 p . As a consequence of (I3.24p . one has7 (n) (£,L) = 
25*( n )(0), where S(' n \k) satisfies the integral equation [k > 0) 



S {n \k) 



g 2 / dtK(V2gk,V2gt) 



cos in 



h 



sinh I 



g 2 / dte-*K(y/2gk,y/2gt)—-rS( n) (t) 
Jo smh § 



(Ins)-" 



(4.8) 



with the symbol |(i ns )-« standing for the component proportional to (lns)~ n . The analysis of this case 
was done in detail in the paper [30]. We report here the main results. 

In order to fulfil condition (I2.10p . the position of the internal holes has to expand in inverse powers 
of Ins: 



u 



(i) 



(4.9) 



n=l 



Introducing the (even) derivatives in zero of the (even) function u(u) = Z'(u), developing them in 
powers of In s, 

-a(u = 0)= ^^(lns)-, (4.10) 

n=-l 



du 2 i 



and imposing the condition f)2.10p for the holes we eventually get the following recursive equation for 
the unknowns a n h, 



p— r+l 



a 



p+l,h 



E- 

.-1 er, 



(-1) 



e n 



[at m ,h 



\.lm 



P-l P-l (I) 



r=1 O'l.-Jp-H-l} m=1 



Ji 



EE^ 



p— r— i+1 



r-1 
(-1) 



e n 



1=0 r=l ^0 {ii,...J p _ r _ i+1 } m=l 



?r(2/i- 1 + L) 



cr. 



(-1) 



p = 



P > 1 
(4.11) 
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where the j m contained in the second term of the right hand side are constrained by the conditions 

p— r+l p— r+1 p—r—l+1 p—i — l+l 

2 im = r + l, Yl m 3m = P+ 1 and the ones in the third term by jm = r , X] m 3m = V~^- 

m=l m=l m=l m=l 

The next step is the Neumann expansion for the functions S^ n '(k) (in the domain k > 0): 

S^(k) = f2 S P(9) Jpi ^ 9k) =► 7 (n) (<^) = V2gS[ n \g) . (4.12) 
P =i 

Straightforward but lengthy calculations, originating from equation (13.261) . lead to the conclusion that 
the Neumann modes Sp l \g) satisfy the syster 



- -(2 P -i) f +x d_t r n ( g ,t)j^(V2gt) _ 2(2p _ 1) ^ ViiA)s<;1(g)i 

2 m=l 

(4.14) 



'0 011111 2 m= l 

In (I4.14p V n (g,t) appears as a coefficient in the high spin expansion 



P(s,g,t) = J2'Pn(gM^sr n (4.15) 



of the internal holes- depending function 



n=l 



L-2 



P(s^,t) = ^[cos^ ) -l] , (4.16) 



h=l 



appearing in (I4.8p . As a consequence of (14. lip . V n (g, t) depends on the various coefficients a m ^ of (14. 9 p 



as 

L-2 



7rr ^ h=1 



^ n?n=l ( a m,h)"" 



^n(»,*) = COSy ^ 1 n n-r+l ■ , > Jm = T ' £ = 71 ' ( 41T ) 

r=l {ji.-Jn-r+l} ^ m=1 Jm ' m=1 m=1 

In order to study the system (14.141) . it is useful to introduce the "reduced coefficients" Sp k \ defined in 
equations (4.23,4.24) of [34] as solutions of the reduced systems 

oo 

Sip\g) = tJ(g)-ApJ2Z2 P ,m(g)(-irS^(g), 

m=l 

(4.18) 

oo 

s { 2 pUg) = tl 1 (g)-2(2p-i)J2z2 P -i, m (9)S i J:\g), 

m=l 



8 We use the notation: 

7 ( s [ + °° dt J n (V2gt)J m (V2gt) 

z n , m {g)= I — ■ ( 4 - 13 ) 
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with the explicit (i.e. not depending on the densities ) forcing terms, 



-OG 



t)=r i dh^Mvm (419) 

2n sinh # 



Indeed, solutions to the systems (I4.14p are linear combinations of the various Sr k \g) with coefficients 
depending on a n ^- In particular, for what concerns 7^(5'), its expression in terms of §[ and a m ^ 
reads 



L-2 

rn—r+l , 



-2tt Si /2) cos — ^ h=1 +1 : , ^ Jm = r , £ mj 



v Z 9 r =l {ji,...,j n _ r+1 } llm=l .7m.- m=l m=l 

(4.20) 

The various a ni /, are written in terms of the densities with the help of (14. lip , leaving eventually ^ n \g) as 
depending on S[ k ^ and a^ q , with r < n — 3. This last property is very important, since it makes possible 
to build up a recursive calculation scheme for the 7^(<?), opening the way to push the computation up 
to the desired order in Ins. 

As an example in paper [3U] we gave the following exact results for the first 7 < - n - ) (g): 

l W (g,L) = 0, (4.21) 



3 



^ 2 \g,L) = V2g - (L - 3)(L - 2)(L - 1) S[%) , (4.22) 



3(a 



r3 JO) 



^(g,L) = -2V2g * \ (L - 3)(L - 2)(L - 1) ffl(y) , (4.23) 



3(a 



vr 2 / a« 3 (a^ 2 



7«C**) = ^{[-^(^-^)(L-3)(L-2)(L-1 



7T <7 



90 (a 



_(L - 3)(L - 2)(L - 1)(5 + 3L(L - 4)) S\ 



(-1)^5 
1 



7T 



^-(L - 3)(L - 2)(L - 1)(5 + 3L(L - 4)) s{ 2) } . (4.24) 



360 (a^\ 

Due to the aforementioned recursive properties, such expressions can be explicitly computed in the 
weak and in the strong coupling limit. Weak coupling expansions are presented in Appendix A of [30J. 
The strong coupling leading term is given in Section 3 of that paper. 

5 High spin and large twist 

We now want to study the anomalous dimension in the limit (II. 7ft : 

s — > 00 , L — )• 00 , j = — fixed . 

In s 
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Since the number of internal holes becomes infinite, we need to treat the sum over them contained in 
(I3.26P by means of the technique discussed in Section 2. Before doing that, we introduce two points ±c 
which separate the internal holes and the Bethe roots: since the number of internal holes equals L — 2, 
the 'separator' c has to satisfy the following relation: 



1 f c 

Z(c) = - / dva{v) = —irj Ins . 

2 J~c 



(5.1) 



Now we can specialise formula (12.191) to the sum 

L-2 



E 

h=l 



COS tu 



CO 
h 



(5.2) 



over the internal holes G [— c, c]. Remember that in [— c, c] no Bethe roots are present: therefore 
the right hand side of (12.191) is zero. We get 



L-2 

£[< 

h=l 



COS tu 



CO 
h 



dv 



(cost-u — l)cr(f ) 



Ti t sin tc 



2ir 6 cr(c) 

7^3^(7(0) sintc + 3t 2 cri(c) cos tc — 3t ^4^r sin tc + toi (c) sinte 



360 



+00 



4^ 



<x(fc) 



(a(c))^ 
sin(t + k)c sin fee 



O 



(Ins) 



k 



7r t sin tc 



6 cr(c) 



77rfV(c) sintc + 3t (Ji(c) costc + tc^c) sintc 



360 



O 



(Ins) 



(5.3) 



where cx m (c) denotes the m-th derivative of the density cr(v) in v = c. This expression has to be inserted 
in (I3.26P and worked out together with condition (I5.ip : we get in general a non-linear integral equation 
for the quantity S(k), related to the Fourier transform of the density of Bethe roots and internal holes 
a(k) through (13^5) . 

This joined analysis of (13.261 15. ip gets simplified in the limit ( II. 7p . Indeed, it follows from the 
structure of (13.261) that, in the case of limit (11.71) . the function S(k) expands as [2U [351 EH] 



s(k) = ( ln *r r E s(r,n) + (( ln 

r=-l n=0 

and, correspondingly, the anomalous dimension behaves as (II. 8p : 

00 00 00 

7 (<7, ^, L) = In s E /nfo)i B + E( ln ^ E ^ + (( ln s )"°°) 



(5.4) 



n=0 



r=0 



n=0 



The functions f n {g), fn(g) are called generalised scaling functions. In particular fo(g) coincides with 
the universal scaling function f(g). 
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Similarly, the separator between internal holes and Bethe roots, c, enjoys the following scaling in 
the limit fTTTj) . 

oo oo 



r=0 n=l 



The constants c^ r,n ^ are connected to the various components cr^a m w hich the density and its deriva- 
tives in zero expand in the limit (U.7p 

^(- = 0)=EE4' n) (l-rT, (5-6) 

r=— 1 n=0 

by means of (15.11) : 

1 f c 

- / dva{v) = —T\j Ins. 

2 J-c 

For instance, using (15.51 15 . 61) in ( 15. ip . we get, for the first c^ r,n ^: 

ct (-i,o)' ^(-i.o)^' - 6 [ -(-i,o)]4 ^(-i^p > 

/ \ rr(°'°) , s ^(0,0) (-1,1) 

~(1>1) ^ ^(1,2) 

C y ' ' = 71- — , : rrr- , C y ' ' = — A1X - 



3 



c (1 ' 3) = 3tt F i— J - -vr 3 — — 1— + 



(-1,0)]4 3 [ a (-l,0)]5 6 [ a (-l,0)]< 



, s [^(0,0)12 

Therefore, once we take the limit (I1.7P and use condition (15. ip . we easily get that equation (I3.26P splits 
in a set of equations, one for every function S^ r ' n >{k). Remember that i5'^ _1 ' )(fc) coincides with the BES 
density (function S^ 1 \k) of last section): therefore, it satisfies the BES linear equation (14.21) . i.e. 

/•+00 j. 

S(- 1 '°)(k) = Ag 2 K{V2gk, 0) - g 2 / dte~* K(V2gk, V2gt)—— T S i - lfi \t) . 

Jo smh 2 

Another particular case is the function S^~ l,1 \k): this satisfies the linear equation 

1 r+00 1 _ p| 

S^(k) = -[1 - J (V2gk)} + g 2 ^ dte"K(V2gk, ^Igt)—^ - 

- g 2 dte-*K(V2gk t V2gt)— —rS^'^it), (5.8) 
Jo smh § 

and determines the generalised scaling function fi(g) = 2S^~ 1,1 \0). Equation (15. 8ft was studied in paper 



The last case that has to be treated separately is r = j = 0. However, in this case, the density is 
obtained from S^(k) ((Ins) contribution at fixed twist, see last section) and S^" 1 ' 1 ^), by means of 
the equality 

5 (0) (fc) = (L — 2)S^ 1 ' l \k) + S^°\k) . (5.9) 
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In the three cases we have discussed up to now, the dynamics of the internal holes is not relevant: for 
the required approximations the internal holes can be supposed all lying at the origin, i.e. the sum (15. 2p 
can be considered as vanishing. 

The particular form of the expansion ( 15. 5p with coefficients given by (15. 7\\ is however of fundamental 
importance for all the remaining cases. They can be studied together, by means of the integral equation 



S<r>")(k) = -i 



dtK(\[2gk, y/2gt) 



cos tu 



(i) 



sinh | 



g 2 / dte-^k(V2gk,V2gt)——rS^ n \t) 
Jo sinh \ 



(Ins)'' 



(5.10) 



in which the sum over the internal holes is evaluated through ( 15. 3 p and the 'separator' c is expanded 
as in (15. 5p . with coefficients depending on the quantities cr^'^ and determined by the use of (15.11) . 

From the methodological point of view, we need to distinguish two cases. The first one covers the 
values r = — l,n > 2 and r = 0, n > 1. In this case only the first term in the right hand side of (15. 3p . 
which is linear in the density, is relevant. Therefore in the final equation non-linearity comes only from 
the non-linear dependence of c on the density and their derivatives in zero: 



S (r ' n) (k) 



-°° dt k(y/2gk,y/2gt) 
sinh | 



dv 
2^ 



(costv — l)a(v) 



(In i 



+ 00 



dte-^k(V2gk, V2gt)—^S^ n \t) . 

sinh 2 



(5-11) 



For instance, using (15 ,7p we get for the first cases: 

0, 



dv 

— (costv — l)a(v 



dv . . , 

— (costv — l)a(v 
2% 



c dv 
c 27 
c dv 
„2ir 



(costv — l)a(v 
(costv — l)a(v 



c dv 

— (costv — l)a(v 

c 2VT 

c dv . 1N . 

— [costv — l)a(v 



2tt 

30 [at" 1 ' )] 5 1 K9) 



4 

7r a 



In s-j 2 



1 

-7T 



t 1 



ttH 4 



dv . , , 

— (costv — l)a(v) 
2tt 



7r 2 [c7 (-l,l)]2 t 2 



\ns-y> 120[(T(- 1 ' )] 4 

.- =0, r = l,2 , 
(o,o )(T (-i,i) „ 



(lns) u 



2[ f7 (- 1 >°)] 4 
dv 
2^ 



7i 



(Ins)" 



(Ins)" 



1 , 4 (0,0) 

1 /7r a\ 
6VT[crM.o)]5 



i 2 tV^ 



30[a(- 1 '°)] 5 
costf — l)cr(w) 



-7T 
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In the remaining cases, i.e. r > 1, the evaluation of the sum over the internal holes (j5.2p involves also 
terms explicitly non-linear in the density. Thanks to formulas (I2.19[ [573]) . however, everything is under 
control and, for instance, for the first values of r, n we get 
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In both cases, next steps are the usual ones and details can be found in [341 [351 EQ]- First, we perform 
a Neumann expansion for the even functions S( r,n \k), in the domain k > 0: 



S^ n \k) = f2s£> n \g) Jp ^ 9k) 
p=i 



(5.14) 



This implies that the generalised scaling functions are expressed as 

fi r \g) = V2gSp n \g) . 



(5.15) 



The Neumann modes Sp r ' n \g) satisfy linear systems and are linear combinations of the "reduced" coef- 
ficients Sp k ^ which are solutions of the systems (14. 18ft . The coefficients driving such linear combinations 
depend on the densities and their derivatives in zero, crljj , with r' < r, n' < n — 1. This property 
allows to find, step by step, exact expressions for the generalised scaling functions fn\g) in terms of 
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a% ' n ^ and S[ K> . For the first of them we get the following results: 
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Explicit expressions in the weak and strong coupling limit for the various fn\g) can be given without 
much ado, because of the iterative structure of the relations which determine them. They can be found 
in [SUES EQ]. 

In particular, strong coupling limit of the anomalous dimension is of interest since it can be checked 
against string theory calculations. For what concerns the function 



(5.16) 



n=0 



performing such a check is not a difficult task, after results by Alday and Maldacena [H]. Introducing 
(at large g) the quantity 

287T 1 _E3_ r „ ( \ 

m(g) = ^-rK\9 4 e ^ 
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l + O 



(5.17) 



in [J4] it was proved that in the limit (JT77J) , when g — > oo, j <C g, with j/m(g) fixed, the quantity 
f(g,j) + j has to coincide with the energy density of the ground state of the 0(6) non-linear sigma 
model with mass gap m(g). When j/m(g) C 1 we are in the nonperturbative regime of the 0(6) non- 
linear sigma model. In this case the energy density can be computed by using Bethe Ansatz related 
techniques. This computation has been systematically performed in [12] • In order to have agreement 
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between our calculations for f(g,j) and computations of [12] for the 0(6) non-linear sigma model, we 
must have that the quantities fl n (g) computed in that paper have to be related to f n (g) by the relation 
f n (g) = 2 n - l Vt n {g). Our results [27] for the strong coupling limit of f 3 (g), f 4 (g), f 5 (gjfl are 



h{g) 
h{g) 
his) 



7T 



, +0 (e~v5 ) 



7T 



12[m(<?)] 2 
vr 4 



5! + 0(1) 



+ 



7T" 



640[m(#)] 3 8[m(g)} 3 



where we used the compact notations 
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and agree with corresponding formulas contained in [12] . 
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6 Conclusions 

In this paper we have discussed our activity on the study of minimal anomalous dimension of twist 
operators in the sl(2) sector of M = 4 SYM at high spin. We preferred to give a general description of 
the methods we used: the reader can refer to the original papers for details and for a complete list of 
our results. 

The main tool we used is integral equation ( 13.261) for the density of Bethe roots. This equation 
coincides with the exact non-linear integral equation equivalent to the ABA equations of the si (2) 
sector if we neglect terms going to zero faster than any inverse power of the logarithm of the spin. 
Equation (I3.26P is linear, apart from the non-linear dependence of the internal holes positions u% on 
the counting function Z(u). These non-linear effects are taken into account by inverting relation (12. 10[) : 
in the high spin limit such inversion is feasible because of the recursive properties of the equations (14. lip 
determining the various u^'. This allows to give explicit expressions (I4.20p for the coefficients of the 
high spin expansion of the anomalous dimension. 

When the twist goes to infinity, the internal holes are described by their density and contribute to 
integral equation (I3.26P introducing explicitly non-linear terms. These non-linear terms are evaluated 
(see Section 5) by using properties and techniques of the non-linear integral equation, which are discussed 
in Section 2. Their contribution to main equation ( 13.261) is obtained by applying formula ( 12.191) and is 
reported in equations (15.121 15.13P . In this case also, recursive properties of the equations determining 
the high spin limit of the density are crucial in order to obtain explicit expressions for the anomalous 
dimension. 



9 In [27] results for the strong coupling limit of f n (g) and checks with 0(6) non-linear sigma model results were 
performed up to n = 8. 
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